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Structural characteristics of positionally-disordered lattices: relation to the first sharp 

diffraction peak in glasses 

J. K. ChristieE S. N. Taraskin, and S. R. Elliott 
Department of Chemistry, Cam,bridge University, Lensfield Road, Cambridge, CBS lEW, UK 

(Dated: February 2, 2008) 

Positional disorder has been introduced into the atomic structure of certain crystalline lattices, 
and the orientationally-averaged structure factor S{k) and pair-correlation function g{r) of these 
disordered lattices have been studied. Analytical expressions for S{k) and g{r) for Gaussian po- 
sitional disorder in 2D and 3D are confirmed with precise numerical simulations. These analytic 
results also have a bearing on the unsolved Gauss circle problem in mathematics. As the positional 
disorder increases, high-fc peaks in S{k) are destroyed first, eventually leaving a single peak, that 
with the lowest-fc value. The pair-correlation function for lattices with such high levels of positional 
disorder exhibits damped oscillations, with a period equal to the separation between the furthest- 
separated (lowest-fc) lattice planes. The last surviving peak in S{k) is, for example for silicon and 
silica, at a wavevector nearly identical to that of the experimentally-observed first sharp diffraction 
peak (FSDP) in the amorphous phases of those materials. Thus, for these amorphous materials at 
least, the FSDP can be regarded as arising from scattering from atomic configurations equivalent 
to the single family of positionally-disordered local Bragg planes having the furthest separation. 

PACS numbers: 61.43.-j 
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Medium-range order (MRO) in the structure of amor- 
phous and liquid systems continues to attract a lot of 
attentionii*^ The nature of the short-range order (SRO) 
in such systems is often accessible to diffraction or 
spectroscopic experiments, the physical interpretation of 
which, in the form of atom-atom pair-distribution func- 
tions, is clear. ^ The more subtle features which contribute 
to MRO, by contrast, are harder both to measure and to 
understand. Of particular interest is the so-called first 
sharp diffraction peak (FSDP), a characteristic feature 
of very many amorphous network materials^'* This peak 
occurs in the structure factor at low values of scattering- 
wavevector transfer, k (= |k|). The real-space informa- 
tion contained in the structure factor can be obtained 
via a Fourier transform, and so the FSDP at low k cor- 
responds to atom-atom correlations at large distances; 
hence the presence of the FSDP is usually thought of as 
a signature of MRO.^ Much effort has been expended 
in trying to deduce the real-space structural features 
which contribute to MRO and the FSDP in various amor- 
phous systemsj^i^iSiL2i2ii£ explanations have been pro- 
posed, for example, based on the packing of voids in 
the material,'''^, the existence of quasi-Bragg planes, ^^ or 
other unknown effectSfS. Gaskell and Wallisifii have stated 
that in many amorphous materials (although not all), 
the wavevector of the FSDP is close to the wavevector 
of a diffraction peak (the lowest) in a related crystalline 
phase, and thus, in real space, to a set of Bragg planes 
(with the largest interplanar spacing). They proposed 
that remnants of these lattice planes exist in the disor- 
dered phase, and that scattering from these planes causes 
the FSDP. 

In this paper, we show that the wavevector of the 
FSDP in certain amorphous materials can be predicted 
from the lattice structure of their crystalline counter- 
parts. Amorphous materials are positionally and topo- 



logically disordered; that is, the atoms have shifted from 
their crystalline positions. However, if the wavevector 
of the amorphous FSDP is close to that of the lowest 
peak in the crystalline structure factor^*' it is reasonable 
to assume that the MRO in the amorphous phase could 
contain structural configurations closely related to the 
crystal. The structure of amorphous materials was mim- 
icked in this work by introducing positional disorder into 
the structure of certain lattices, and the behaviour of 
the structure factor and atomic-density pair-correlation 
function was investigated analytically and numerically. 

Positional disorder was introduced into the crystal lat- 
tices by shifting each atom independently a random dis- 
tance away from its crystalline position. In each of the 
X-, y- and ^-directions, the shifts were randomly chosen 
from a Gaussian distribution centred on zero, indepen- 
dently of the shifts in the other directions. The half width 
a of the Gaussian was thus a measure of the amount of 
positional disorder introduced. No restriction was placed 
on the closeness of approach of atoms (the results were 
not qualitatively different with a suitable cut-off). 

If a is small, the resulting disordered structure is rather 
like the original crystal; if a is large, the atomic posi- 
tions are located effectively randomly and little memory 
remains of the underlying lattice. Figure ^ shows the 
structure of the diamond lattice (Fig. ^) as an increas- 
ing amount of positional disorder is introduced. The un- 
derlying lattice is still very clear at cr = O.lri (Fig.^), 
where ri is the crystalline nearest-neighbour distance. 
For a — 0.2ri (Fig.QJ;), some remnant of the {111} layers 
is still visible. Only faint vestiges of the crystal structure 
persist at cr = 0.3ri (Fig.[Tli). 

Gaussian positional disorder, as studied here, has a 
physical analogue: it corresponds to thermal motion of 
atoms about their crystalline positionsiii In this work, 
the displacement of the atoms from their crystalline po- 
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FIG. 1: (Colour online) Gaussian positional disorder in the di- 
amond lattice with values of disorder parameter a: (a) 0, (b) 
O.lri, (c) 0.2ri, (d) 0.3ri, where r\ is the nearest-neighbour 
distance. All models have the same orientation. Atoms in 
three consecutive {111} lattice planes perpendicular to the 
page are darkened to show the gradual destruction of the lat- 
tice planes with increasing positional disorder. 



sitions can far exceed that which would occur thermally 
in a real crystal, but the mathematical formalism is iden- 
tical. In particular, the effect of thermal motion on the 
structure factor ^(k) = 1 + ;^ X^i ^j^i exp(ik.rij), where 
N is the number of atoms, and r^j = Yj — r^ is the 
position vector for the separation between atoms i and 
j, is well-knowniiSiH The crystalline Bragg peaks have 
a weight proportional to iVji4 thermal disorder reduces 
their weights by the Debye- Waller factor, e~ " (with- 
out broadening the peaks) and a low-amplitude diffuse 
scattering intensity occurs. In terms of the (orientation- 
ally averaged) crystalline structure factor S'(fc)'^'"^**, the 
disordered structure factor for Gaussian positional dis- 
order is given by>i2iH S{k) = 1 + e-''^''\S{kyy''* - 1). 
The Debye- Waller factor destroys the high-A; peaks in the 
crystalline structure factor fastest as positional disorder 
increases (Fig. |21). It is clear that, as a increases, there 
will be a range of disorder in which only one peak makes 
a substantial contribution to S{k), the rest being negli- 
gible; this peak will be the crystalline peak at the lowest 
value of k. It is knowniS that the FSDP often has a 
very similar wavevector to that of the lowest-fc peak in 
^^^^crj/st. ]^gj.g .^g ggg ^Yia,t this peak is the most resis- 
tant to positional disorder. The lowest-A: peak in S{k) 
corresponds to the set of Bragg planes which are fur- 
thest separated in real space; they are also the set of 
planes with the highest atomic density. Figure |21 shows 
the calculated S{k) for positionally-disordered diamond 
lattices, with N — 97336 atoms, for various values of a. 
The peaks decrease in height in very close accord with 
the theoretical Debye- Waller decrease. For a > 0.3ri, 
only the {111} peak at kri = 4.71 makes a substantial 
contribution. 

The structure factor S{k) is related by a Fourier trans- 
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FIG. 2: (Colour online) Orientationally-averaged structure 
factor per atom, for Gaussian-positionally-disordered dia- 
mond lattices (with 97336 atoms) with different values of dis- 
order parameter, o, where ri is the nearest-neighbour crys- 
talline distance. Only the lowest-fc peak, at fcri = 4.71, con- 
tributes substantially for a > 0.3ri. 



form to the pair-correlation function g{r) (the orienta- 
tional average of 5 (r) = (iVpo)"^ Ej Ej^* (^ (^ " fy), 
where po is the average density). Thus, the effect of 
positional disorder is to convolve the (5-function peaks 
in the crystalline ^(r) with the Fourier transform of the 
Debye- Waller factor. For Gaussian disorder, this gives: 
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where d is the dimensionality of the system. 

In principle, g{r) can be evaluated by summing eq. 
n over all crystalline interatomic distances for a given 
lattice, and orientationally averaging the result. For the 
2D square lattice with crystalline lattice positions r^ — 
rix'x. + HyY, g(r) in eq. ^can be orientationally averaged 
to give: 
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where Io{z) is a modified Bessel function of the first kind. 
For the 3D simple cubic lattice, with crystalline lattice 
positions r^ = n^jX + Uyy + n^z, g{r) is given by: 
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where n — \ n^ + n't + ^21 ^^^ ^hc sums in eqs. I2l3l ex- 
clude the case where all indices ni are zero simultane- 
ously. It has not been possible to evaluate these sums 
analytically. In retrospect, this failure is unsurprising. 
More than 150 years ago, Gauss first posed what is known 
as the Gauss circle problemji^ given a regular square lat- 
tice of points, and a circle of radius r centred on one of 
the points, how many lattice points lie within the cir- 
cle? The exact answer remains unknown, despite much 
attention. ^^ An exact expression for g[r) with a width a 
of zero would solve this problem. The equivalent problem 
is also unsolved in more than two dimensions. 

Since it has not been possible to obtain g[r) analyti- 
cally, it was calculated numerically by simulating large 
(^ 10^ atoms) crystalline models of 2D square and 
3D simple cubic (SC), face-centred cubic (FCC), body- 
centred cubic (BCC), diamond, and cristobalite (Si02) 
structures, which were then Gaussian positionally disor- 
dered with various values of a. Periodic boundary con- 
ditions (PBCs) were used, and it was necessary to use 
large models in order to increase the number of allowed 
k-vectors, particularly for low k. PBCs also imply that 
g{r) is not reliable for distances beyond half of the (cubic) 
box size L; data are not shown for values of r > L/2. 

For the crystal (with cr = 0), g[r) consists of a se- 
ries of (5-functions at the various interatomic distances. 
For small a (for a < 0.2ri - not shown here), the S- 
functions broaden into Gaussiangil and merge, giving a 
complicated form. However, as a is further increased, 
g{r) takes on a simple oscillatory form with a single fre- 
quency, having an envelope that decays spatially as a 
power law (see Fig.|2}. The 3D lattices were all found to 
have g{r) decaying as r~^ , whilst the 2D square lattice 
was found to have g{r) decaying as r^^/^. The reason 
for this is clear: the 3D Gaussian g{r) (eq. 13) has an ex- 
plicit r~^ dependence, and the asymptotic behaviour of 
the modified Bessel function in the 2D g{r) (eq. [2Jl has 
an r~^^^ dependence. The simulated 5(r)'s (Fig.lSJ agree 
very closely with the theoretical expressions given by eqs. 
I2I3I If g{r) consists of oscillations with a single signifi- 
cant frequency of period D, then there is only one signif- 
icant peak in the structure factor S{k), at a wavevector 
of fci = 2Tr/D. A heuristic argument for why only a sin- 
gle period is observed in the oscillatory behaviour of g{r) 
at sufficiently large a draws on the Rayleigh criterion for 
resolving overlapping broad peaks, such as those that oc- 
cur in the sums of Gaussians in eqs. I2I3I two peaks are 
just resolvable when the full width of the peaks is com- 
parable to their separation. Thus, for sufficiently large 
a, only contributions from lattice planes with the largest 
separation combine to give resolvable oscillations in g{r). 

The periods of the oscillations, D, were extracted by 
assuming the empirical form g{r) = 1 + Ar~^ sin{kir + (f)) 
for the 3D lattices, and g{r) = 1 + Ar^^^^ sm{kir + (f>) 
for the square lattice, and performing a least-squares 
minimisation with respect to the variables A, ki and (p. 
Non-oscillatory data for r < 2ri were excluded from the 
fit. Averages were then taken over twenty realisations of 




FIG. 3: (a) g{r) for a Gaussian-positionally-disordered di- 
amond lattice witli disorder parameter a = 0.3ri (solid), 
and the best-fit curve (dashed), (b) The same data plotted 
double-logarithmically with the best-fit r~^ decay envelope. 
(c) g{r) for a Gaussian-positionally-disordered square lattice 
with a = 0.25ri (solid), and the best-fit curve (dashed), (d) 
The same data plotted double-logarithmically with the best- 
fit r~^" decay envelope. Both best-fit curves exclude the 
non-oscillatory data at low r. 



each disordered model. For all the lattices studied, D is 
equal to the spacing between the furthest-separated lat- 
tice planes (see Tabled). We conclude that these single- 
frequency oscillations in g{r) correspond to the Fourier 
transform of the single remaining peak in S{k). 

We now show that positionally-disordered structures 
of appropriate crystal lattices can predict the position 
of the FSDP in various amorphous systems. Real amor- 
phous materials show several peaks in S'(fc), all of which 
comprise contributions from many interatomic distances. 
The corresponding g(r) has several aperiodic peaks, and 
the single-period oscillation shown in Fig. O would mod- 
ulate such peaks, giving rise to the FSDP. Recently, 
"time-frequency" wavelet-transform analysis has been 
applied to the study of disordered structures to elucidate 
the interatomic-distance contributions to different peaks, 
particularly the FSDPJU&ia We show below that the 
FSDP in S{k) for some amorphous materials could be 
due to scattering from structural configurations equiva- 
lent to the last surviving family of local, disordered Bragg 
planes in a positionally-disordered structure. 

Silicon crystallises in the diamond structure, and its 
amorphous form has been well studied as a canonical 
example of a topologically-disordered solidi^ It also ex- 
hibits a strong FSDP at fci = 2.00A-i.Sa From Table E 
the period of the last surviving extended-range oscilla- 
tions in g{r) for the positionally-disordered diamond lat- 
tice is the distance between successive {111} planes, i.e. 
the apex-basal-plane distance in an SiSi4 tetrahedron. 
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Lattice 


Furthest planes 


Interplanar distance 


Period of oscillations, D 


SC 


[100] 


ri 


(1.0000 ± 0.0003) n for a = 0.25ri 


FCC 


[111] 


yfri = 0.8165ri 


(0.8167 ± 0.0001) n for a = 0.25ri 


BCC 


[110] 


Tfri = 0.8165ri 


(0.8165 ± 0.0001) ri for a = 0.25ri 


Diamond 


[111] 


fn = 1.3333ri 


(1.3333 ± 0.0001) n for a = 0.35ri 


Square 


[10] 


ri 


(0.99993 ± 0.00001) n for a = O.Sn 



TABLE L Families of lattice planes with the largest interplanar distance for certain crystal structures; ri is the nearest- 
neighbour distance in each case. For the values of disorder parameter, a, given, g{r) consists of oscillations having only a single 
frequency, the periods of which, D, were found from a best-fit procedure, as described in the text. 



spending to an oscillation with a period oi D = 3. ISA, 
and hence fci — 27r/3.13A= 2.0lA~'^, in excellent agree- 
ment vi^ith the experimental value aboven^S 

The FSDP is often observed in AX2-type glassesji*^ 
the most well-known example being vitreous silica, v- 
Si02. A crystalline counterpart of v-Si02 is the tetrag- 
onal a-cristobalite structure. ^^ Real v-Si02 is known^ 
to consist of fairly ordered Si04 tetrahedra, linked in 
a corner-sharing manner in a disordered network; we 
made no attempt to keep the tetrahedra ordered, and 
Gaussian disordered the atomic positions of cristobalite 
without restriction. Hence, these models are not a par- 
ticularly good representation of real v-Si02. Neverthe- 
less, single-frequency extended-range oscillations in g(r) 
were found in positionally-disordered cristobalite struc- 
tures for 0.35 < cr/ri < 0.6. The period of this oscilla- 
tion is D ~ 4.O7A, corresponding to a peak in S{k) at 
ki ~ 1.54:A~^, in very good agreement with the experi- 
mental value of the FSDP wavevector at fci = 1.52A^^.^^ 

In conclusion, Gaussian positional disorder has been 
introduced into a variety of crystalline lattices, specifi- 
cally, the 3D diamond, FCC, BCC, SC and the 2D square 



lattices, as well as the cristobalite crystalline polymorph 
of silica. As this positional disorder increases, the high-/c 
Bragg peaks in the structure factor are destroyed first, 
eventually leaving only the lowest-fc peak in S{k). Ex- 
amination of the orientationally-averaged real-space pair- 
correlation function of such disordered structures shows 
power-law-damped, single-frequency, extended-range os- 
cillations, with a period, D, equal to the spacing between 
the most distant (most areally dense) lattice planes. This 
period is equal to D — 27r/fci, where fci is the wavevector 
of the first peak in the structure factor. The atomic posi- 
tions of crystalline counterparts of silicon and silica have 
been disordered. The wavevector of the FSDP in real 
amorphous silicon and in vitreous silica, is nearly iden- 
tical to that of the lowest-A; peak in the structure factor 
of the positionally-disordered crystalline counterparts, 
and hence at sufficiently large distances from an arbi- 
trary origin atom, the atomic structure of a topologically- 
disordered amorphous material appears to be similar to 
that of the corresponding positionally-disordered crystal. 
JKC would like to thank the Engineering and Physical 
Sciences Research Council for financial support. 
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